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© The invention is comprised of a method which 
allows the user to convert several data layers from 
one GIS into another. The user creates a template 
map of the area of interest. The template is an "all 
points" map which, when graphed, would appear as 
a black polygon identical in shape to the area of 
interest. Using the first GIS, the template is con- 
verted to a file in bit-mapped format. All zero bits in 
this file are identified and their position in the file 
noted. This information is referred to as the trans- 
form function. For each data layer to be converted, 
the first GIS is used to display a map of the proper 
extent. This display is converted to a file in bit- 
mapped format and the transform function used to 
delete non-data bits which correspond to the zero 
bits identified in the template file. The resulting, 
modified, bit-map is imported into the second GIS. 
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FJELD OF INVENTION 

The present invention relates to computer aid- 
ed geographic information systems (GIS). More 
particularly, this invention relates to an improved 
method of converting data from one GIS to another, 
maintaining a consistent system of georeferencing. 

BACKGROUND OF THE INVENTION 

A geographic information system is an informa- 
tion system designed to work with data referenced 
by spatial or geographic coordinates. For example, 
on a very basic level, a contour map of an area ov 
land (that is, a map in which contour lines indicate 
the elevations at specific geographic points) couid 
be considered a GIS. If, to the contour map, we 
added a street map, a sewer map, and an aerial 
photograph, then we would have a GIS composed 
of four data sets: each set being referred to as a 
data layer or data plane. 

In an automated GIS, the data layers will be 
stored in a data base system and a set of tools will 
be included to enter, manipulate and analyze the 
data, and then display the result either on a screen 
or as hardcopy printout. 

To create an automated GIS, it is necessary to 
first identify and gather data. Typical sources of 
such data are municipal maps, government survey 
maps, aerial photographs and other publicly avail- 
able data sources. The data must be extracted 
from these sources and then be manipulated so 
that it may be entered into the GIS. In addition, the 
data must be georegistered or georeferenced. That 
is, the spatial data must be referenced to a coordi- 
nate system such as Universal Transverse Merca- 
tor, State Plane, or Latitude/Longitude. 

The geo-registration process is mathematically 
as follows: 

x = f!(X,Y) 
y = f 2 (X,Y) 

where 

(x,y) = distorted coordinates in some co- 
ordinate system, 
(X,Y) = correct coordinates in a selected 

reference grid 
fi ,f2 = transformation functions 
In its simplest form, where the distorted coordi- 
nates are taken from a scale drawing and the 
selected reference grid is a flat surface marked by 
latitude and longitude, the transformation functions 
are linear equations of the type: x = AX + BY + 
C, where A, B and C are constants. In the case of 
aerial photographs, where the intent is to represent 
the curved earth surface on a flat map, and where 
the angle of the photograph, as well as the pho- 



tographs representation of a curved surface in a 
two dimensional plane must be taken into account, 
the transformations are considerably more complex 
but are well known ir. the an. 
5 The standard method by which data is 

georegistered against selected reference coordi- 
nates (i.e. of identifying the required transformation 
functions) can be described as follows: 

1. select a point from the spatial data available 
10 (distorted coordinate system) and determine the 

equivalent point in the selected reference grid; 

2. express both selected points in terms of pairs 
of coordinates and compute the transformation 
function which reflects the relationship between 

15 the pairs; 

3. choose another point and repeat steps 1 and 
2. 

Repeating the process for all points in the 
spatial daia set would be prohibitively costly in 

20 both time and effort. Accordingly, the process is 
repeated only for a selected sample of points. The 
method by which the points are chosen (the sam- 
pling methodology) can be taken from a variety of 
well-known techniques such as nearest neighbor, 

25 bilinear, or cubic convolution. 

It is not uncommon for the "true" transforma- 
tion function to vary over the area of interest. For 
example, in an aerial photograph, the curvature of 
the earth results in a different mapping of points to 

30 the two dimensional picture. Accordingly, the 
"true" transformation function varies depending on 
the distance from the observation point to the pic- 
tured surface. As another example, consider a dia- 
gram of water mains. If the purpose of the diagram 

35 was to indicate the points of interconnection, then 
the representation of the mains between such inter- 
connection points would not necessarily be con- 
sistant in the sense that there would likely not be a 
constant ratio of mapped lines to the physical dis- 

40 tance between points. The "true" transformation 
function tc georeference such a diagram would 
consist of a collection of linear functions, none of 
which were necessarily related to any other. 

Where the "true" transformation varies over 

45 different parts of a map, the transformation cal- 
culated from only a subset of points may have 
certain inaccuracies. This wili result in distortions 
when the calculated transformation function is used 
to georeference the entire set of spacial data 

so points. 

Various schemes for automating the process of 
georeferencing exist in the art. Most of the tech- 
niques focus on reducing the distortion by choos- 
ing a large number of coordinate pairs (i.e. rubber- 

55 sheeting techniques.) The current state of the art is 
described in the book Remote Sensing and Image 
interpretation , Lillesand and Kiefer, Second edition, 
1987, Chapter 10, Section 10.2. 
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Often GIS users require information embodied 
in data in more than one GIS. Since each GIS 
typically keeps its data in a proprietary format, 
there are very few commercially available conver- 
sion programs which allow a user to import data 
from one or more source GISs to a single target 
GIS. 

In the prior art, conversion from one GIS to 
another was done by taking the physical output of 
each of the required data layers that resided in one 
GIS and mapping that output onto the second GIS. 
To describe this process more specifically, the, 
following notation will be used: Data Layer (a,b) will 
refer to a bth layer of GISa. Accordingly, assume 
that a user wishes to use Data Layer (1,1), Data 
Layer (1,2), and Data Layer (1,3), all of which are 
found in GIS1, together with Data Layer (2,1) which 
is found in GIS2. Further assume that GIS2 con- 
tains a particularly useful data manipulation tool so 
that it is mandatory that the GIS1 Data Layers 
(1,n) be converted to GIS2, rather than the more 
simple step of converting Data Layer (2,1) to GIS1 
format. 

In this situation, the user. would take the spatial 
information of Data Layer (1,1) (typically by using 
the tools of GIS1 to create a graphic output such 
as a physical map) and georegistering this informa- 
tion against Data Layer (2,1), using the same tech- 
niques that would have been used had Data Layer 
(1,1) and Data Layer (2,1) been a surveyor map 
and a desired coordinate system, respectively. (In 
other words, using the techniques described in the 
preceding paragraphs.) This process would be re- 
peated for Data Layer (1 ,2) and Data Layer (1 ,3). 

Since each of the Data Layers (1 ,n) contained 
different data, (different points and lines would be 
represented on the map associated with each Lay- 
er), the set of points the user selected from Data 
Layer (1,1) to register against points from Data 
Layer (2,1) would necessarily be different from the 
set chosen from Data Layer (1,2) to be registered 
against Data Layer (2,1). This means that the time- 
consuming process of repeatedly choosing differ- 
ent coordinate pairs must be executed for each 
layer. Further, as previously discussed, the choice 
of a particular subset of coordinate pairs can intro- 
duce certain distortions in the georegistration. 
Since the sets of coordinates were different for 
each Data Layer (1 ,n), the distortions in the 
georegistration would be different for each layer. 
This result is that the imported Layers would be 
registered imperfectly and would not overlay prop- 
erly. 

OBJECTIVES OF THE INVENTION 

Accordingly, it is an objective of the invention 
to provide a means of converting spatial data from 



one GIS to another GIS so that the registration 
between the layers of the first GIS is not lost and 
the problem, of improper overlay is minimized. 

It is a further objective of the invention to 
5 provide a method of converting multiple layers of 
spatial data from one geographic information sys- 
tem to another whereby the problem of determining 
multiple transformation functions is avoided, thus 
providing consistent distortions in the transformed 
10 . layers, said consistent distortions being more easi- 
ly compensated for. 

These and other objects* features and advan- 
tages of the invention are achieved by the tech- 
nique described herein. 

75 

SUMMARY OF THE INVENTION 

The invention is comprised of a method which 
allows the user to convert several data layers from 

20 one GIS into another. The user creates a template 
map of the area of interest. The template is an "all 
points" map which, when graphed, would appear 
as a black polygon identical in shape to the area of 
interest. Using the first GIS, the template is con- 

25 verted to a file in bit-mapped format. All zero bits in 
this file are identified and their position in the file 
noted. This information is referred to as the trans- 
form function. For each data layer to be converted, 
the first GIS is used to display a map of the proper 

30 extent. This display is converted to a file in bit- 
mapped format and the transform function-is used 
to delete non-data bits which correspond to the 
zero bits identified in the template file. The result- 
ing, modified, bit-map is imported into the .second 

35 GIS. . ^ • * 

■ , -.* . t - . . .* * ■ 
BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 
FIG. 2 



FIG. 3 



is a flow diagram of the invention, 
is a graphic representation of a se- 
lected geographic extent, the tem- 
plate for such extent, and the bit-map 
associated with such template, 
is pseudo code used in creating the 
transformation function. 



DETAILED DESCRIPTION OF THE INVENTION, 

The invention described herein, namely, a 
50 method for georegistration of imported bit-mapped 
spatial data into a Geographic Information System 
(GIS), finds application in any computer-based GIS. 
In its preferred embodiment, the invention was ex- 
ecuted on a host data processing system com- 
55 prised of a workstation, an IBM RISC System/6000 
Model 530 with 64 Megabytes of memory and a 
1.6 Megabyte hard-drive, having connected to it a 
display screen IBM 6091-19 and keyboard and 
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mouse. The RISC System/6000 used the AIX Ver- 
sion 3 operating system under which were ex- 
ecuted tv/o GIS software programs. Examples of 
such GIS programs are ARC/INFO™, a program 
provided by ESRI (Environmental Systems Re- 
search Institute) and GRASS™, a GIS program 
provided by United States Construction Engineer- 
ing Research Laboratory. 

The invention will be described with reference 
to the problem of taking information from GIS1 and 
mapping it into GIS2. In the preferred embodiment, 
the information consisted of that part of three data 
layers which was applicable to a fixed geographic 
area or extent. For example, the three data layers 
could be respectively land use, water, and a con- 
tour map, and the geographic extent or area could 
be a 20 block area of Denver, Colorado. 

FIG.1 is a flow diagram of the invention. The 
steps herein described may be embodied in a 
program stored in the RISC System/6000 work- 
station, 

The method cf the invention begins with the 
selection of the geographic extent of interest, Siep 
10 of FIG. 1. In step 100, a template is created to 
correspond to this geographic extent. The template 
is a two-dimensional representation of the entire 
geographic extent of interest to the user. Thus, a 
template is a polygon whose shape is congruent to 
the area of interest. More specifically, both the 
geographic extent of interest and the template are 
areas bounded by line segments, each of which 
are connected to another at its end points. The 
length of the line segments of the geographic ex- 
tent are in a constant ratio to the length of - the 
corresponding line segments of the template and 
the angle between adjacent line segments is the 
same for both the geographic extent and the tem- 
plate. Referring to FIG. 2, it will be seen that in the 
case under discussion, the template would be a 
rectangle (20), the length of whose sides bear a 
mathematical relationship to the distance of the 
area cf Denver that is under analysis (10), that is, 
the area which extends from longitude L1 to L2 and 
latitude L3 to L4. 

Further, the template is an "all points" tem- 
plate. That is, each point within the template is 
considered to be a data point, those points outside 
of the template area are not. Referring to FIG. 2, it 
will be seen that the one embodiment under dis- 
cussion, the template (20) created is a two dimen- 
sional mask, the interior (21) of which is black 
(represented by the shading in FIG. 2), the exterior 
(30) of which is white. Thus, for this situation 
"data" 's equivalent to "black" and "non data" to 
"white." 

Step 110 is the determination of the coordinate 
pairs of the template extent, that is, the coordinate 
pairs of the end points of each of the line segments 
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bounding the geographic extent. In the example 
under discussion, this would consist of four data 
points, corresponding to the four corners (22, 23, 
24, 25) of the template. In FIG. 2, the coordinate 

5 pairs in question would be (L1,L3), (L2, L3), (L1.L4) 
and (L2,L4) respectively. 

Step 120 comprises displaying the Template 
by using the display functions of the GIS. Exam- 
ples of such functions are the command "arcs" in 

io ARC/INFO and "d.rast" in GRASS. 

In Step 130, the information cn the screen is 
captured in a file format. A number of techniques 
are available to execute this step. Many GISs pro- 
vide a screen print option (since most users of a 

15 GIS require hard copy of the maps they create 
from the system). For example, in ARC/INFO, un- 
der Arcplot the use of the "screensave" option 
saves the current screen display as a raster image. 
There are other programs, unrelated to GIS, which 

20 allow the user to save the screen area in a file. For 
example, the X-Window Dump command "xwd" 
will save screen areas. The output of this step is 
what is referred to as a bit-map file, FILE T, which 
represents the template as displayed. As is known 

25 in the art, a bit map is a representation of char- 
acters or graphics by individual pixels or points of 
light, dark or color, arranged in row (horizontal) and 
column (vertical) order. Each pixel is represented 
by either one bit (simple black and white) or up to 

30 32 bits (for high definition color). In the embodi- 
ment under discussion, black is represented by a 
"1" bit and white by "0." A representational picture 
of the bitmap for template (20) is shown in FIG. 2, 
. (50). 

35 Step 140 is the identification in FILE T all bits 

in the bit-map which are not "1" (or "data"). In this 
step, all records which are not bit-map data, for 
example, header records or header fields, and all 
parts of records which are not bit-map data are 

40 identified. The location of each of the "0"s or non- 
bit-map data is thus identified. 

FIG. 3 represents pseudo code for a program 
which isolates the data from non-data in a template 
(300). As can be seen from FIG. 3, the code 

45 processes the bit-map first by row (31 G) and then 
by column (350). 

The output of Step 140 is the transform in- 
formation which will be saved and applied to all 
subsequent layers. The geographic template extent 

so (in the example depicted in FIG. 2, ((L1,L3), 
(L2.L3), (Li,L4) and (L2.L4)) will also be saved at 
this point. 

Step 150 consists of using GIS1 to display a 
map which consists of Data Layer (1,1) information 
55 for the template extent. 

Step 160 consists of creating a bit-map file, 
FILE L1 , in the same way that FILE T was created 
from the Template in Step 130. 

4 
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Step 170 consists of using the Transform Tem- 
plate against FILE Li. That is, the transform in- 
formation is used to compress FILE L1 by dropping 
all the bits in the positions which correspond to the 
zero bits in Step 140. The output of this step is the 
transformed file, FILE L1\ 

Step 180 consists of converting FILE LV to 
ASCII data using programs such as GRASS'S 
"r.in.ascii" command and the template extent iden- 
tified in Step 110. 

As can be seer> from FIG. 2, Step 190, Steps 
150 through 180 are repeated for Data Layer (1,2) 
and Data Layer (1 ,3), creating files File L2 1 and File 
L3\ All three transformed files (in this example, 
FILE L1\ FILE L2' and FILE L3') are imported to 
GIS2. 

As will be clear to those skilled in the art, an 
important step in the invention and improvement 
over the prior art is the step of recognizing that a 
template map can be created and exploited for all 
succeeding coverages or layers to be transferred. 
Since the transform function is created only once, 
for all points, by creating the function for the tem- 
plate map and then applying that function to all 
subsequent map layers, all subsequent map layers 
can be imported without re-selecting coordinate 
pairs, or recomputing transformation functions. This 
results in perfectly co-registered map layers in the 
second GIS, as well as correct geo-regi strati on to a 
coordinate system. In addition, if a distortion is 
introduced in the creation of the transformation 
function, the distortion is carried systematically 
throughout each of the data layer transforms. Sys- 
tematic errors to those skilled in the art can either 
be easily corrected or, for some applications, even 
ignored. 

As discussed, FIG. 3 is an example of a portion 
of the instructions which can provide the function 
corresponding to Step 140 of FIG.1, namely, to 
convert the bit-map of the template to a transform 
function. For the embodiment under discussion, a 
program with the function of the pseudo code 
shown in FIG. 3 was written in the C programming 
language as described in the publication "The C 
Programming Language" by Brian W. Kernighan 
and Dennis M. Ritchie, 1978, and also "X Window 
Systems Programming and Applications with Xt w 
by Douglas A. Young, 1989. Programs with this 
function can be written in any other conventional 
program language, as, for example, Fortran. 

It is within the scope of the invention, that the 
order of the steps of the instructions in FIG. 3 can 
be altered, or other steps added, without changing 
the fundamental nature of the invention. It is simi- 
larly within the scope of the invention, that the 
order of the steps of FIG. 1 can be altered, or other 
steps added (as for example, providing for error 
routines), without changing the fundamental nature 



of the invention. 

In as much as this invention outlines a method 
for georegistering the bit-mapped display output of 
one GSS to be provided as input to another GIS, the 

5 technique described herein solves several impor- 
tant user requirements. For example, in the ab- 
sence of a common data exchange format, this 
invention gives users of GIS products the ability to 
preserve a significant investment in their data (typi- 

io cally around 80% of total system implementation 
costs).' 

The invention has been described in reference 
to specific embodiments. Although a specific em- 
bodiment of the invention has been disclosed, it 

75 will be understood by those having skill in the art 
that other embodiments, variations and modifica- 
tions to the herein described specific embodiment 
can be made without departing from the spirit and 
scope of the invention. For example, it will be clear 

20 to those skilled in the art that the invention may 
also be used to map one or more data layers from 
GIS1 and one or more data layers from GIS2 into 
GIS3. Accordingly, it is not intended that this inven- 
tion should be limited except as indicated by the 

25 accompanying claims. 

Claims r 

1. In a graphic computer system including a dis- 
30 play and a plural-ity of geographic information 

systems (GIS), each of said GIS containing a 
plurality of data layers, each data layer com- 
prised of georeferenced data, a method for 
transferring georegistering data between a plu- 
35 rality of GIS while preserving georeferencing 

and maximizing accuracy; said method com- 
prising the steps of: 

a; selecting a geographic area, said area's 
boundary comprised of a plurality of line 
40 segments, each of said line segment being 

adjacent to two other of said line segments, 
and each line segment having two end 
^ points, each of said end point being shared 
with one of said adjacent line segments. 
45 Each of said geographic area end points 

having a coordinate pair associated with it, 
said coordinate pair identifying the geo- 
graphic location of said end point, 
b. creating a template, comprising a two- 
50 dimensional polygon whose boundary is 

comprised of a plurality of line segments, 
each of said line segment being adjacent to 
two other of said line segments, and each 
line segment having two end points, each of 
55 said end point being shared with one of said 

adjacent line segments. Said template being 
congruent in shape and area to said geo- 
graphic area and each of said template end 
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points being associated with one of said 
geographic area end points; 

c. determining for each of said template end 
points, the coordinate pair of said asso- 
ciated geographic area end point; 

d. using a first GIS to display said template 
on said display screen, each interior point of 
said template on said display being of a first 
data type and each exterior point of said 
template being of a second data type; 

e. creating from said template display, a 
template file, said template file comprised of 
a plurality of bits, each of said bits identified 
with a specific point on said screen; 

f. creating a transform function, said trans- 
form function consisting of identifying in 
said template file all bits which are of said 
second data type and recording the location 
of such bits in said file; 

g. using said first GIS to displaying a por- 
tion of a first data layer from GIS1 on said 
display screen as a first map, said portion 
comprised of data georeferenced to said 
selected geographic area; 

h. creating from said displayed first map a 
first file, said first file comprised of a plural- 
ity of bits, each of said bits identified with a 
specific point on said screen; 

i. applying said transform function tc said 
first file to create a first transformed file, 
said application comprised of identifying in 
said first file all bit positions corresponding 
to positions with 0 bits in said template file 
and deleting said identified bits; 

j. converting said first transformed file to the 
input forma* of a second GIS; 
k. using said first GIS to displaying a portion 
of a second data layer from GIS1 on said 
display screen as a second map, said por- 
. tion comprised of data georeferenced to 
said selected geographic area; 
I. creating from said displayed second map 
a second file, said second file comprised of 
a plurality of bits, each of said bits identified 
with a specific point on said screen; 
m. applying said transform function to said 
second file to create a second transformed 
file, said application comprised of identify- 
ing in said second file all bit positions cor- 
responding to positions with 0 bits in said 
template file and deleting said identified 
bits; and 

n. converting said second transformed file 
to the input format of said second GIS. 

A method as in claim 1 wherein, following step 
n of Claim 1, said first transformed file is 
displayed overlaying said second transformed 
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file are displayed using said second GIS. 

3. A method as in claim 1 wherein said first data 
type is a "1 "-bit and said second data type is 
a "0"-bit. 

4. A method as in claim 1 wherein said coordi- 
nate pairs comprises the latitude and longitude 
of said end points. 

5. A method as in claim 1 wherein said display 
on said screen is comprised of points of light, 
dark or color, said points being referred to as 
pixels, said points being arranged in rows and 
columns. 

6. A method as in ciaim 5, wherein said first file 
created in step h is a bit-map, said bit map 
comprised of a collection of bits, each bit 
associated with one of said pixel of said dis- 
play, said value of said bit determined by 
whether said pixel is a point of light, dark or 
color, and said position of said bit in said bit- 
map being determined by said row and col- 
umn of said associated pixel. 

7. A method as in claim 6, wherein said second 
file created in step 1 is a bit-map. 

8. A method as in claim 1 wherein said second 
data layer is a data layer of a third GIS. 

9. A method as in claim 6 wherein said second 
data layer is a data layer of a third GIS. 

10. A method as in ulaim 7 wherein said second 
data layer is a data layer of a third GIS. 
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FIG. 1 
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FORMAT, CREATING A BIT MAP FILE Ln 
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CONVERT FILE Ln* TO DATA FORMAT ACCEPTABLE TO GIS2 
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RETURN TO STEP 145 



\ 

BNSDOCID: <EP 0619S54A2_I_> 



j=P^S19 554 A2 



FIG. 2 



L2 

1 SPACIAL 
DATA 




0000000000000000000000000000000000 BIT-MAP 

OOO00000OOOOO000OO00OOO000OOOOCK30O 

0000000000000000000000000000000000 



000000001111 
000000001111 
000000001111 
000000001111 
000000001111 
000000001111 
000000001111 



111100000000000 
111100000000000 
111100000000000 
111100000000000 
111100000000000 
111100000000000 
111100000000000 



50 



0000000000000000000000000000000000 

oooooooooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooooooo 



BNSOOCID: <EP 0619554A2J_> 




Eui opaisches Patentamt 
European Patent Office 
Office europeen des brevets 



l lUllilliyili.il 111 I!,! 

© Publication number: 0 619 554 A3 



EUROPEAN PATENT APPLICATION 



© Application number: 94102228.7 
(§) Date of filing: 14.02.94 



© int. CI. 6 : G06F 15/72 



® 


Priority: 17.03.93 US 32211 




Inventor: Berger, Blaine H. 






2941 Bowline Place 


® 


Date of publication of application: 




Longmont, CO 80501 (US) 




12.10.94 Bulletin 94/41 




Inventor: Kressin, Mark S. 








7598 Extate Circle 


® 


Designated Contracting States: 




Niwot, CO 80503 (US) 




DE FR GB 




Inventor: Smith, Bret P. 








742 Cameron Lane 


® 


Date of deferred publication of the search report: 




Longmont, CO 80501 (US) 




05.07.95 Bulletin 95/27 












Representative: Teufel, Fritz, Dipl.-Phys. 




Applicant: International Business Machines 




IBM Deutschland Informationssysteme 




Corporation 




GmbH, 




Old Orchard Road 




Patentwesen und Urheberrecht 




Armonk, N.Y. 10504 (US) 




D-70548 Stuttgart (DE) 



© A method for geo-registration of imported bit-mapped spatial data. 



CO 
< 

m 
in 

o> 

CD 



CL 
UJ 



© The invention is comprised of a method which 
allows the user to convert several data layers from 
one GIS into another. The user creates a template 
map of the area of interest. The template is an "all 
points" map which, when graphed, would appear as 
a black polygon identical in shape to the area of 
interest. Using the first GIS, the template is con- 
verted to a file in bit-mapped format. All zero bits in 
this file are identified and their position in the file 
noted. This information is referred to as the trans- 
form function. For each data layer to be converted, 
the first GIS is used to display a map of the proper 
extent. This display is converted to a file in bit- 
mapped format and the transform function used to 
delete non-data bits which correspond to the zero 
bits identified in the template file. The resulting, 
modified, bit-map is imported into the second GIS. 
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FIG. 3 



Read Template Map File 
Determine number of rows (NROWS) 
Determine number of columns (NCOLS) 
ZerovariablesR1,R2,C1,C2 
Set Boolean flag R_FLAG TO FALSE 
310 For (row=0; row < NROWS; row ++) 

Read one row of data 

If R_FLAG = FALSE and entire row is zero then R1=row 

Else R_FLAG = TRUE 
If R_FLAG = TRUE and entire row is zero then 
DO 
R2 = row 
row = NROWS 
END 

END 

350 Read one row of data at row = R1 
Set Boolean flag C.FLAG to FALSE 
For (column = 0; col < NCOLS; column++) 

If C_FLAG = FALSE and column is zero then C1 = column 

Else CFLAG = TRUE 
If C_FLAG = TRUE and column is zero then 
DO 
C2 = column 

column = NCOLS C1 C2 
END 
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